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Abstract 

We propose a spin foam model of four-dimensional quantum gravity which 
is based on the integration of the tetrads in the path integral for the Palatini 
action of General Relativity. In the Euclidian gravity case we show that the 
model can be understood as a modification of the Barrett-Crane spin foam 
model. Fermionic matter can be coupled by using the path integral with 
sources for the tetrads and the spin connection, and the corresponding state 
sum is based on a spin foam where both the edges and the faces are colored 
independently with the irreducible representations of the spacetime rotations 
group. 
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1 Introduction 



The approach of defining a quantum theory of gravity by using a path integral 
quantization has been revitalized by the appearance of the idea of spin foams 
[1]. A spin foam model can be described as a lattice gauge theory for a BF 
theory, and although a BF theory is a topological theory, the Palatini action 
of General Relativity (GR) can be represented as a constrained BF theory, 
where the two-form B is a wedge product of the spacetime tetrade one-forms. 
This then leads to the idea that the GR path integral could be defined as 
a modification of the path integral for a topological theory. This was the 
approach used for the construction of the Barret-Crane (BC) models [2, 3], 
which culminated when a finite partition function for GR was constructed for 
any non-degenerate triangulation of the spacetime manifold [4]. However, it 
was soon realized that one can obtain several finite BC models with different 
convergence properties [5]. This ambiguity is a problem because it is still 
not clear which one of these has GR as the classical limit. The source of 
the ambiguity is the fact that the edge amplitudes of the dual two-complex 
cannot be fixed in the BC quantization procedure, which then leads to many 
possible models. 

Another problem with the BC type models is that it is difficult to couple 
matter, especially fermions, because matter fields couple to the tetrades and 
it is often impossible to rewrite the matter actions coupled to gravity as 
functionals of the B and matter fields only. Several proposals for coupling 
matter have been made [6, 7], but none of these proposals are based on 
manipulations of the matter plus gravity path integral, so that it is not clear 
what is their physical relevance. 

These problems of the BC approach suggest that one should try to find a 
spin foam model which is based on the integration of the tetrade fields in the 
GR path integral. Such an approach should be feasible because the Palatini 
action is quadratic in the tetrads, so that the path integral over the tetrads is 
Gaussian. In this paper we show that the effective gauge theory obtained by 
the integration of the tetrads can be defined as a spin foam model. In section 
two we describe the construction of the path integral for pure gravity, while 
in section three we describe the construction of the path integral for gravity 
plus matter, including the cosmological constant term. In section four we 
present our conclusions. 
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2 Pure gravity path integral 



We consider first the case of pure gravity, and we start from the Palatini 
action 

S=[ e abcd e a Ae b AR cd = [ (e 2 R)d A x , (1) 

where M is the spacetime manifold, e a are the tetrad one-forms, R ab = 
du ab + uj a c A uj cb is the curvature two- form, uj ab is the spin connection one- form 
and e a bcd is the totally antisymmetric symbol (eom — !)■ The corresponding 
path integral can be rewritten formally as 

Z = J VuVee^M {e2R)dix = J Vuj (det R)' 1/2 , (2) 

where (deti?)~ 1//2 denotes the result of the integration of the tetrads. 

The formal expression (2) suggests that one may try to define Z on a 
triangulation of M as 

z = I n^n(detF/)- i/2 = / n^n a (^/) . (3) 
if if 

where A x = f t u, g t = e A \ detF = (e abcd F ab F cd ) 2 , 

9f = e F f= II 9i , A(g f ) = (det F f y^ , (4) 

l&df 

and the indices / and / stand for the edges and the faces of the dual two- 
complex of the triangulation. The group function A(g) should be gauge 
invariant, so that we take 

A( 9 ) = J:A(A)xa(9) , (5) 

A 

where Xa(s') is the character for an irreducible representation (irrep) A, and 
the sum is over all irreps of a given category (finite-dimensional or unitary). 
This then implies that 

A(A) = f G dgxK{g)A{g) . (6) 
By using the formula 

IdgD^ig) ■ ..D^(g) = E« W (c££ W )* , (7) 
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where C^Zaf ^ are the components of the intertwiners t for the tensor prod- 
uct of four irreps are the corresponding representation matrices, 
we will obtain a state sum of the form 

z=En A ( A /)n4,(A/,i,) , (8) 

A/A / v 

where the vertex amplitude A v is given by the evaluation of the pentagon 
spin network, which in the SU (2) case is known as the 15j symbol. The state 
sum (8) is of the same form as in the case of the topological theory given by 
the BF action; however, the weights we put on the faces are not dim A/ but 
the functions A (A/). 

These new weights are given by the integrals which are generically diver- 
gent, due to det-F/ = configurations, so that some kind of regularization 
must be used. In order to do this, let us write the Lie algebra element Ff as 

Ff — Ef • K + Bf • J , (9) 

— * — * 

where K are the boost generators, while J are the spatial rotations genera- 
tors. The so(4) Lie algebra is a direct sum of two so(3) algebras, and a basis 
of this decomposition is given by 

J± = l(J±ZK) , (10) 

where £ = 1 in the Euclidian case and £ = % in the Minkowski case. 
From (10) it follows that 

det F={E-Bf = ^{{E + f-{E_f) , (11) 

— * — * 1 — * —* 

where E± = B ± jE. Note that E± are real in the Euclidian case, while in 
the Minkowski case are complex conjugates (E + = EJ). The group function 
A(g/) is then given by the expression 

AG/) = |((£+) 2 - (i-) 2 )" 1 , g = e?+-& +S -*-=g + g- . (12) 

Let E± = \E±\, then the weights A(Aj) can be calculated from (6) and 
(12). Since D^ A \g) is a direct product of two SU (2) representations D^+\g + ) 
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and D^~\g_), we obtain 

A ( A ) _ f d 3 E + sin 2 (E+/2) f d 3 E_ sin 2 (£_/2) 



rs 8tt 2 El 7r3 8tt 2 £ 2 
sin(j+ + \)E + sin(j_ + 4 



sin(£ + /2) sin(£_/2) £(£ 2 - El ' 
We have used the SU (2) group integration measure 



(13) 



where g = e E ' a is a spin | representation matrix [8], and the trace formula 

X {9) ~ ME/2) ■ U5j 

Since d 3 E = E 2 sin 9 dE d6 dip, and there is no angular dependence in the 
integral (13), we can put d 3 E = 4nE 2 dE, so that 



£vr 2 A(A) = /_ dE + » + (E + ) /_ ^/'V > ( 16 ) 



(17) 



00 dE_fi_(E_) 

o — rr " rv - r/ y -E 2 — -E 2 

where fJ,±(E) = sin(£'/2) sin(j± + |)E. The integral 

/■+« dE_n_{E_) 
Jo E\-E 2 _ 

is divergent, and has to be regularized. By using 

sin(£/2) sin (j + ^ E = cos(jE) - cos(j + l)E , (18) 
we obtain I — Ij — Ij+i, where 

dE cos{kE) 



I, 

Since 



k = Jo E\-E 2 ■ (19) 

1 dEe ikE . . 

h= 2 Re Lw^ ' (20) 
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it has a form of a one-dimensional quantum field theory (QFT) propagator. 
One can then regularize I by the QFT methods, i.e. by replacing by 
+ ie (e > 0) and by performing a contour integral along the real axis and 
a semi-circle in the upper half-plane [9]. After taking the e — * limit, one 
obtains 

sm(kE + ) 



h = -*^ET • (21) 

so that 

/■oo dE 

-ZnA(j,l)=J o — [cos(jE)-cos(j + l)E}[smlE-sm(l + l)E] .(22) 

In order to evaluate (22) we will study the integral 

^fe~ zE , zeC, e>0. (23) 

This integral can be rewritten as 

I(z) = [%e~ zE (24) 
00 dx 

— e~ zex = E^ze) . (25) 

Ob 

Since 

E 1 (z) = - 1 -]ogz-Y, ' , , \argz\<7i , (26) 

where 7 is the Euler-Mascheroni constant [10], then 

I(z) = — In e — 7 — In \z\ — i arg z + 0(e) , |arg,2|<7r . (27) 



Let 



then 



/■oo HE 

I(jJ)= J o -^cos(jE)sin(/E) , (28) 



/OO WR' 

— [sin(j + /)£ - sin(j - Z)£] (29) 

= Ilim/m[E 1 (^(j + /) e )-E 1 ( ? (/-j>)] (30) 

= \* + \arg(i(l-j)) = ^e(l-j) . (31) 



Therefore 



A{j,l) = ~mi-j)-6(l-j + l)-e(l-j-l)] . (32) 

This implies that the non-zero coeficients are the ones with I — j = or 
I — j = ±1, whose values are given by 

A(j,J-l) = ! , A0\j) = m^ , A (J , J + 1) = ^1 , (33) 

where # — $(0). The natural values for 0(0) are and 1, although 9(0) — \ 
gives a more symmetrical A(j,l). In any case, one obtains a weight which 
is concentrated around the simple irreps so that the model can be 

considered as a generalization of the Barrett-Crane model. 



3 Path integral with sources 

Including matter and the cosmological constant term requires the evaluation 
of the path integral 

Z = J VeVujV-iPexp^i J M ((e 2 R) + X(e 4 })d 4 x + iS m [-ilj,e,uj}^ , (34) 

where A is the cosmological constant, S m = f M d 4 x£ m and C m is a function 
of the tetrads, spin connection, matter fields ip and their derivatives. Note 
that in the case of spin-half fermions C m is a polynomial in e and u given by 

S m = J m e abcd e a A e b A e c A ^ (7" (d + \uj rsl r l s ) + mej ^ , (35) 

where 7 a are the Dirac gamma matrices and m is the fermion mass. Hence 
the path integral (34) can be evaluated at least perturbetively. 

In the case of scalar and vector bosons C m is not a polynomial function 
of the tetrads. However, in the vector boson case, when the Yang-Mills 
action is put on the simplicial complex, it is possible to obtain a polynomial 
approximation suitable for the spin foam evaluation [11]. 

When C m is a polynomial of the fields and their derivatives, the path 
integral (34) can be evaluated perturbatively via 

Z= lim e ^/M<W«) 4 >^s m [-«/«x,-«/«,-^ Zo ^ >J>x ] , ( 36 ) 
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where 

Z [J,j, X ] = J VeVuV^e l hM e2R )+ J ^ ab +^ a +^) di - (37) 

is the generating functional. Z is essentially the gravitational path integral 
with the sources since the matter integration in (37) gives a delta function 
S(x)- The tetrade path integral is Gaussian, so that we need to define 

Z [J,j] = J ^e^M^^^det^-^^-^M^-JM^^^-)/?- 1 ^,^'^)^ 

(38) 

This expression can be defined on a triangulation of M along the lines of 
the J = j = case. However, when the sources are present a more intricate 
state sum will appear. Guided by the expression (38) we will define Z as 

Z (J,j)= [U d 9iK9i,Ji)Il A (9f,Je) , (39) 
i f 

where 

l*(g h Jt) = e^toV , A(g f , Je ) = A(g f )e-^ F i 1 ^ 4 , (40) 

and the subscripts e and e denote two edges of the triangle dual to the face 
/■ 

The function /j, is not a gauge invariant, and therefore it will be expanded 

as 

^g h J l ) = Y,KKJ l )D^ ) {g l ) . (41) 
Hence the coefficients //(A, J) are matrices given by 

/i(A,J)= / dgD\g)^g,J) . (42) 

J G 

The function A will be expanded as a gauge invariant function, so that 

A = E A ( A / ? ^)XA / (^/) • (43) 

A/ 

Hence 

A(A,j)= / dg X A(g)A(g,j) . (44) 
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The group integrations in (39) can be performed by using the analog of 
the formula (7) for the tensor product of five irreps. One then obtains a state 
sum 

zoUj)= e n A (A/,i e )n^ A ^ j on^( A /' A ^o • ( 45 ) 

A/,Ai,t; / I v 

This is a novel spin foam state sum, because it involves a dual 2-complex 
whose edges and faces are independently colored with the irreps of the group. 
Zo(J,j) can be understood as an amplitude for a Faynman diagram given 
by a five-valent graph whose edges carry the irreps A; and the loops carry 
the irreps A/. The edges have propagators /x(A/, J{) and each loop carries a 
weight A(Aj, j t ), while the vertex amplitudes A v are given by the evaluation 
of the pentagon spin network with five external edges, where the internal 
edges carry Aj irreps, while the external edges carry A; irreps. 

The propagators /i(A, J) will be given by the integrals (42), where fi(g, J) = 
e i(a + -j++a--j-) ^ g = e a+-j++a--j- and A = (k, I). In order to see what kind of 
integrals are these and what kind of regularization should be used, we will 
evaluate their traces, which are given by 

Trn{K,J)= I dgTrfi(g,J)xA(g) ■ (46) 

J G 

We then obtain 

<-°° duj 



Tr/i(k,l,J±) = / — ± I2 + (lo + ) / de + sme + e^+ J + cos0 + 

JO 4:71 J 

/ ^ /U _(cu_) / ^_sin^e^- J - cos9 - . (47) 

JO 47T JO 



This integral can be easily evaluated by using the formulas from the previous 
section, so that we obtain 

TrMk,l,J ± ) = 0[J + ~ k) ~ e(J + - k ~ l) B(J ' - S - 0[J - . 

j _)_ j — 

(48) 

Note that this function is non-zero only for the spins belonging to the 
intervals [J± — 1, J±], so that the sum over the A; irreps is a finite sum. 
However, the weights (48) are not analytic functions in the vicinity of J = 0, 
so that in the limit J — > the partition function will diverge. In order 
to regularize this, we will calculate (47) perturbatively in J by expanding 
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e iu>jcos6 - m Baylor series around J = 0. The /i will be then a power series 
in J, with the coefficients given by the integrals 

jj 1 /*00 I f'7T 

Cn (k) = — dusm(u/2)sm( y k + -)u}( y uJ) n dd sin 9 cos n 9 . (49) 
n! Jo 2 jo 

After performing the 9 integration, we obtain 

Cn(k) = p— / dw[cos(kuj) — cos(k + l)u]uj n , (50) 

where n is an even number. This integral can be regularized by performing 
an analytic continuation of 

roo 

/ e~ auJ u n du = n\ a- 71 - 1 , Rea>0 , (51) 
Jo 

to Re a = 0. This gives 

c 2n {k) = ^[{k + I)"**" 1 - A;- 2 - 1 ] , (52) 

so that 



n+ 2 



TrfjL(k,l,J ± )= J2 c 2 n(k)c 2m (l)Jl m J- n ■ (53) 

m,n>0 

The weights A(A, j) are the character expansion coefficients of the group 
function 

" Vdet/ j ' (M) 

where we have used F' 1 = -§F*(det F)' 1 ' 2 . By using the formulas from 
the previous section, we obtain 

A(g,j) = (El - E 2 _y l exp [i ^+ ' |+ ± f 2 ^ ) , (55) 

where Q± are quadratic functions of j and j' . The weights A will be then 
given by the integral 

00 dE 4 



^±/i + (£ + ) Jd9 + sm9, 



47T 

cos g + +Q_£'- cos ( 



r°° dE r exp u b^-^ / 
i ~4t~^~ Jo Jg-V Z - (56) 



This integral can be evaluated perturbatively in Q±, which is sufficient 
for our purposes, since at the end of calculation of the partition function the 
sources j have to be set to zero. One can use the Taylor series 

( .Q+E+ cos6> + + Q-E- cos6L\ 
exp % 



El - E 2 



co -n 



= W=W- j ' (5?) 

and because of the angular integrations only the even powers of n contribute. 
The corresponding integrals can be regularized by the same procedure as in 
the previous section, i.e. by replacing (E+ — El)^ 1 with (E+ — E 2 _ + ie) -1 , 
e > 0, integrating along the contour in the upper half-plane and then taking 
the e — > limit. 
We then obtain 

A(A,j) = A(A)+a + (A)Q 2 + + a_(A)Q 2 _ 

+ P + {A)Q% + {3_{A)Qt + i{A)Q 2 + Ql + 0{Q & ) , (58) 

where A (A) is given by the equation (32). For example 

dE 4 



, . . I dE + r 



d6+ sin 6 4 



Jo 4?r v ' Jo {El - El + ief v ' 

It is easy to see that a + (k, I) = —«_(/, k) and one can show that a-(k, Z) is 
given by a convergent integral 

a_{k,l) = -j-^j a dEsm{E/2)sm{{k + l/2)E)[f l {E)-f l+l {E)l (60) 

where 

/K^) = -^cos(/E) + ^^sin(/E) . (61) 

Note that in the expansion (58), as well as in the expansion (53), only 
the even powers of the ± components appear, which means that the 50(4) 
invariance is preserved. More explicitely 

J± = - J ab Jab ± £ abCd Jab-Jed , Q± = -^f^fab ± t"**** fab fed , (62) 
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where f a b = \{jaj\, — jbj' a )- Since f a b is associated to a triangle, there will be 
three source vectors ji, j'2 and j'3 for each edge of the triangle, and therefore 
we will take a symmetric expression 



4 Conclusions 

We have formulated an approach to obtain a spin foam model of quantum 
gravity based on the integration of the tetrade fields in the path integral for 
the Palatini action for General Relativity. We have determined the weights 
for the spin foam state sum in the Euclidian gravity case, and one can see 
that the tetrade model is a generalization of the BC model since all the irreps 
of 5*0(4) are included in the state sum but the weights are non-zero only for 
the irreps close to the simple irreps. Beside allowing for the coupling of the 
fermionic matter, the tetrade model is an improvement over the BC model 
since all the simplex weights can be determined (in the BC case the edge 
amplitudes are arbitrary). 

Matter and cosmological constant can be coupled by introducing the 
sources for the tetrades and the spin connection, which results in a new 
type of the spin foam state sum. It is a sum of the amplitudes for the space- 
time two-complex where the edges and the faces carry independent labels. 
If one thinks about a spin foam as a world-sheet of a spin network, then 
including sources corresponds to putting additional spins on the verticies of 
the spin network. This is close to the picture where matter corresponds to 
the external edges of an open spin network, an approach which was used in 
[6] to couple matter. However, the formalism of [6] requires open spin foams, 
while for the tetrade model only closed spin foams are used. 

In order to have a physical model, the simplex weights should be defined 
for the Minkowski case. One way to do this would be to perform an analytic 
continuation of the Euclidian weights (32) such that £ — > i. In this approach 
one would work with the same category of representations as in the Euclidian 
case, i.e. finite-dimensional SU (2, C) x SU(2, C) representations, so that 
the Minkowski weights will be the Euclidian weights times the appropriate 
factors of %. An alternative approach would be to use the category of unitary 
SL(2, C) representations. In this case the irreps are infinite-dimensional and 
can be labeled as (j,p) where 2j e Z + and p e R+. The corresponding 




(63) 
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weights will be determined by the integrals (6), (42) and (44). 

An important next step is to study the convergence of the state sum in the 
Euclidian and the Minkowski case. Note that if the Euclidian state sum turns 
out to be divergent, it can be regularized by passing to the quantum group 
at a root of unity, which is usually done in the case of topological spin foam 
models. However, since our model is non-topological, using a quantum group 
regularization is not necessary, so that one can use alternative regularizations, 
for example a gauge fixing procedure for spin foams [12]. 

As far as the semiclassical limit is concerned, this is still an unsolved prob- 
lem for quantum gravity spin foam models. The difficulty is that in the case 
of non-topological models the partition function state sum is triangulation 
dependent. This means that it is not a diffeomorphism invariant, as it should 
be. Formally one can sum over all the triangulations in order to obtain a 
triangulation independent expression; however, it is still not known how to 
do this. Perhaps adapting the techniques of the dynamical triangulations 
approach [13] could help with this problem. 

An alternative approach to the problem of finding the semiclassical limit 
would be to find a smooth manifold limit. This would require studying 
triangulations with very large number of simplexes, and hopefully one could 
extract an effective diffeomorphism invariant action. However, a technique 
must be developed in order to do this, perhaps something analogous to the 
methods used for the study of the two-dimensional Ising model near the 
critical point. 

The simplest thing one can do in order to better understand the model 
is to analyze the state sum for the simplest triangulation of the four-sphere 
given by the hexagon graph 2 . Clearly, more work is necessary in order to 
resolve these issues, but the model is promising because it was obtained from 
the discretized path integral for the Palatini action. 
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